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Abstract 

Chromosomal inversions disrupt recombination in heterozygotes by both reducing crossing-over within inverted regions and 
increasing it elsewhere in the genome. The reduction of recombination in inverted regions facilitates the maintenance of 
hybridizing species, as outlined by various models of chromosomal speciation. We present a comprehensive comparison of 
the effects of inversions on recombination rates and on nucleotide divergence. Within an inversion differentiating Drosophila 
pseudoobscura and Drosophila persimilis, we detected one double recombinant among 9,739 progeny from Fi hybrids 
screened, consistent with published double-crossover frequencies observed within species. Despite similar rates of exchange 
within and between species, we found no sequence-based evidence of ongoing gene exchange between species within this 
inversion, but significant exchange was inferred within species. We also observed greater differentiation at regions near 
inversion breakpoints between species versus within species. Moreover, we observed strong "interchromosomal effect" 
(higher recombination in inversion heterozygotes between species) with up to 9-fold higher recombination rates along 
collinear segments of chromosome two in hybrids. Further, we observed that regions most susceptible to changes in 
recombination rates corresponded to regions with lower recombination rates in homokaryotypes. Finally, we showed that 
interspecies nucleotide divergence is lower in regions with greater increases in recombination rate, potentially resulting from 
greater interspecies exchange. Overall, we have identified several similarities and differences between inversions segregating 
within versus between species in their effects on recombination and divergence. We conclude that these differences are 
most likely due to lower frequency of heterokaryotypes and to fitness consequences from the accumulation of various 
incompatibilities between species. Additionally, we have identified possible effects of inversions on interspecies gene 
exchange that had not been considered previously. 
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Introduction 

In the study of speciation, a crucial question is how species 
persist despite gene flow. One proposed solution is that 
chromosomal inversions partition the genome into regions 
protected from gene flow by reducing recombination over 
long stretches (Noor et al. 2001; Rieseberg 2001). 
Although reduced recombination in inversion heterozy- 
gotes has been documented repeatedly within species 
(Roberts 1976; Ishii and Charlesworth 1977), few studies 
estimate recombination rates in interspecies hybrids het- 
erozygous for inversions. Further, by focusing on the 
importance of inversions in the reduction of recombination 
within inverted regions, recent work has largely overlooked 
the major global increases in recombination rate in inver- 
sion heterozygotes (but see Portin and Rantanen 1990), 



known as the "interchromosomal effect" (Schultz and 
Redfield 1951). Rather than determining the effects of 
the interchromosomal effect on global nucleotide variabil- 
ity, emphasis has instead been placed on identifying the 
mechanism of this process (Joyce and McKim 2010). 

We examine here the extent to which inversions differen- 
tiating species affect recombination rates 1) within inverted 
regions, 2) at inversion boundaries, and 3) throughout the 
remainder of the genome. We focus our research on hybrids 
between Drosophila pseudoobscura and Drosophila persi- 
milis, a classical model system for studying chromosomal 
inversions, hybridization, and speciation (Dobzhansky 
1 937). Drosophila pseudoobscura is found across Western 
North America extending north to Canada, east to Central 
Texas, and south to Central America. Drosophila persimilis 
maintains a smaller range restricted mainly to the western 
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US and is contained within the range of D. pseudoobscura 
(Dobzhansky and Sturtevant 1938; Dobzhansky and Epling 
1944). These species are morphologically identical; how- 
ever, male hybrids are sterile (Dobzhansky 1936). They 
diverged approximately 0.5-0.85 Ma (Aquadro et al. 
1991; Hey and Nielsen 2004; Leman et al. 2005); however, 
parts of the genome carry a signature of more recent 
hybridization (Machado et al. 2002; Hey and Nielsen 
2004; Machado et al. 2007; Kulathinal et al. 2009). Despite 
this ongoing gene flow, these taxa were diagnosed as dif- 
ferent species with three major paracentric inversions largely 
distinguishing them — two on the X chromosome and one 
on the second chromosome (Dobzhansky and Epling 1 944). 
The D. pseudoobscura arrangement of the right arm of the 
X (hereafter XR) is also present in populations of D. persimilis 
exhibiting meiotic drive with skewed sex ratios, referred to 
as Sex Ratio (SR) D. persimilis (Sturtevant and Dobzhansky 
1936). Additionally, both species segregate multiple 
arrangements among individuals on chromosome 3 
(Dobzhansky and Sturtevant 1938; Powell 1992). 

To examine how inversions affects recombination rates 
within inverted regions, we focused on the largest (XR) inver- 
sion because we expect higher potential for double-crossover 
events and gene flow between the two arrangements 
(Navarro et al. 1 997). We evaluated the double-crossover rate 
in female hybrids between D. pseudoobscura and D. persimi- 
lis across this —12.5 Mb inversion through a direct assay. We 
further compared divergence of different karyotypes be- 
tween-species and within-species for the XR-chromosome 
inversion using published genome sequences of D. pseu- 
doobscura (Richards et al. 2005) and D. persimilis (Clark 
et al. 2007), and an assembled genome sequence we gener- 
ated for D. persimilis possessing the SR arrangement. 

We also examined the effects of inversion heterozygosity 
on recombination rates elsewhere in the genome. To this 
end, we quantified recombination at the boundaries of in- 
versions by estimating the amount of recombination be- 
tween markers outside, but near, the breakpoints of each 
inversion. These regions have higher rates of divergence, 
similar to markers inside inverted regions (Machado et al. 
2007; Noor et al. 2007; Kulathinal et al. 2009; McGaugh 
et al. forthcoming), suggesting that they experience less re- 
combination than neighboring regions. Finally, we calculate 
the effect inversions have on enhancing recombination rates 
throughout the rest of the genome by comparing recombi- 
nation rates in hybrids to a published recombination map of 
D. persimilis (Stevison and Noor 2010). 

Expected Rates of Exchange within Inverted 
Regions Based on Within-Species Inversion 
Polymorphisms 

The major prerequisite of modern chromosomal speciation 
models (for recent review, see Faria and Navarro 2010) is 



that crossover products are rarely or not recovered between 
inverted segments. In Drosophila, single crossovers between 
heterokaryotypes of a paracentric inversion result in nonvi- 
able gametic products, which are subsequently shunted to 
the polar bodies prior to oviposition (Sturtevant and Beadle 
1 936; Carson 1 946; Coyne et al. 1 993; Navarro et al. 1 997). 
Therefore, any exchange between heterokaryotype regions 
is likely due to either gene conversion or double-crossover 
events. Due to the strong effect of interference acting over 
distances as long as 8-10 Mb in Drosophila (Weinstein 
1918; Foss et al. 1993; Fitzpatrick et al. 2009; Stevison 
and Noor 201 0), only inversions with large recombinational 
distances (>20 cM) are expected to achieve an observable 
rate of double-crossover (Navarro et al. 1997). 

Empirical analysis of exchange in inversion heterozygotes 
has followed two approaches: controlled crosses to measure 
the rate of double-crossover and sequence analysis of tar- 
geted regions within an inversion to infer historical exchange. 
Empirical estimates between segregating inversions in Dro- 
sophila have observed an average of 10~ 4 double crossovers 
within a single generation using phenotypic mutants located 
near the center and ends of inversions spanning 30-80 cM 
(Levine 1956; Ishii and Charlesworth 1977). This double- 
crossover rate is lower than recent empirical studies using mi- 
crosatellite markers within inversions spanning 7.7-1 1.7 Mb 
along the O chromosome in D. subobscura, which found one 
double-crossover among 391 individuals in one inversion het- 
erokaryotype cross (~10~ 3 ) and none in two other crosses of 
inversion heterokaryotypes (Pegueroles et al. 2010). In com- 
parison, the double-crossover rate across similar regions in ho- 
mokaryotypes is approximately 5-10% (Gruneberg 1935; 
Spurway and Philip 1952; Novitski and Braver 1954; Robbins 
1974; Pegueroles etal. 2010), showing that double-recombi- 
nation products are not produced as frequently as expected in 
inversion heterozygotes. Sequencing approaches have en- 
hanced our understanding of differentiation within inverted 
regions by being able to detect double-crossover events or 
historical gene flux due to conversion within an inverted re- 
gion (Rozas and Aguade 1 994; Wesley and Eanes 1 994; Has- 
son and Eanes 1 996; Betran et al. 1 997; Laayouni et al. 2003; 
Schaeffer and Anderson 2005; Nobrega et al. 2008; White 
et al. 2009). These studies confirm predictions by theoretical 
models for flux rates to be highest near the center of inverted 
regions (Navarro et al. 1997; Kirkpatrick and Barton 2006; 
Feder and Nosil 2009). 

Despiteour understanding of the ratesof double-crossovers 
in inversions differentiating populations of Drosophila, we 
do not have a strong understanding of how inversions that 
differentiate species might differ in rates of exchange from 
inversions segregating within species. There are three reasons 
toexpect exchange rates might differ in interspecies hybrids: 1 ) 
inversion heterozygotes segregating within populations are 
more abundantthan those between species, which are depen- 
dent on the rate of premating isolation, 2) mechanistically, 
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levelsof sequencedivergence influence rates of crossover, pre- 
dicting lower levels of exchange when sequence divergence is 
higher between homologous chromosomes (Modrich and 
Lahue 1 996), and 3) inversion heterozygotes between species 
will likely suffer some fitness consequences because of alleles 
conferring differential adaptation or incompatibilities within 
them, allowing selection to influence the detectable rate of his- 
torical exchange (Noor et al. 2001; Kirkpatrick and Barton 
2006). Although this latter difference may not influence the 
rate observed in a single-generation cross, it may impact the 
observed exchange rate in a sequence-based analysis. For 
this reason, we used direct cross analysis to calculate dou- 
ble-crossover rate in a single generation despite a lack of ev- 
idenceof ongoing exchange in a recentsequence-based test of 
migration within the inverted region of XR (Stevison 201 1). 
Furthermore, by comparing both SR D. persimilis and 
D. pseudoobscura sequences (both carrying the ancestral 
XR arrangement) to the D. persimilis sequence, which carries 
the derived XR arrangement, we examine how exchange 
within the XR inversion varies in intraspecific versus interspe- 
cific comparisons. 

Recombination Suppression of Inversions Extends 
beyond Breakpoints 

Markers outside inverted regions, but near the breakpoints, 
tend to have heightened divergence, similar to markers 
found within inversions (Machado et al. 2007; Noor et al. 
2007; Kulathinal et al. 2009). This effect is consistent with 
the observation of persistent recombination suppression in 
inversion heterozygotes outside the inversion, near bound- 
ary regions (hereafter inversion boundaries) (Dobzhansky 
and Epling 1948; Spurway and Philip 1952; Maynard Smith 
J and Maynard Smith S 1954; Ortiz-Barrientos et al. 2006; 
Kulathinal et al. 2009; Pegueroles et al. 201 0). This pattern is 
not unique to Drosophila (White and Morley 1955; Roberts 
1976; Stump et al. 2007; Strasburg et al. 2009), and similar 
effects of extended divergence at centromeric boundaries 
suggest that this pattern is not unique to inversion bound- 
aries (Carneiro et al. 2010). Our systematic analysis of 
recombination rate at inversion boundaries aims to refine 
the previous estimates from Kulathinal et al. (2009) for 
the extent of recombination suppression at inversion bound- 
aries in hybrids between D. pseudoobscura and D. persimilis. 

Inversions as Global Recombination Modifiers 

In addition to their proposed effects on species maintenance, 
chromosomal inversions increase recombination rates signif- 
icantly throughout the rest of the genome (Sturtevant 1 91 9; 
Schultzand Redfield 1951). Inversions alter crossover rates on 
the same chromosome (intrachromosomal effect) and on 
other chromosomes (interchromosomal effect). Due to mul- 
tiple inversions in this system, we cannot tease apart intra- 
chromosomal effect and interchromosomal effect in this 



study, so for simplicity, we refer to them collectively hereafter 
as ICE. Despite renewed interest in inversions and their impor- 
tance in speciation (see recent reviews: Hoffmann and Riese- 
berg 2008; Brown and O'Neill 201 0; Faria and Navarro 201 0; 
Jackson 2011), advances in fine-scale recombination map- 
ping technology, and the observation of extensive variation 
in recombination rate from humans to yeast (Stephan and 
Langley 1998; Gerton et al. 2000; Cirulli et al. 2007; Coop 
and Przeworski 2007; Kulathinal et al. 2008; Rockman and 
Kruglyak 2009; Stevison and Noor 2010), there has been rel- 
atively little research done on the role of ICE in contributing to 
nucleotide variation. Here, we map recombination rates on 
the second chromosome in single-generation hybrids (hetero- 
zygous for three inversions, see Introduction) and compare 
these measures with the fine-scale homokaryotype recombi- 
nation map on the same chromosome recently published in 
D. persimilis (Stevison and Noor 2010). 

Materials and Methods 

Single-Generation Estimate of Inversion Crossover 
Rate in Hybrids 

Genome lines of D. pseudoobscura, Mesa-Verde, CO 2-25 
(MV2-25, San Diego stock number #14011-0121.94) and 
D. persimilis, Mount Saint Helena 3 (MSH3, San Diego stock 
number #1 401 1 -01 1 1 .49), were crossed to generate hetero- 
zygous F, hybrid females. These females were backcrossed to 
MV2-25 males to generate approximately 1 0,000 progeny to 
screen for crossovers along the XR inversion. 

Using microsatellite or indel markers at the breakpoints 
and center of the XR inversion (spanning 12.5 Mb), we 
tested for double crossovers in Ft hybrids between D. pseu- 
doobscura and D. persimilis. The rate of gene conversion is 
often assumed in recombination models to be uniform 
across an inversion, whereas the rate of double-crossover 
is expected to be highest near the center, due to the role of 
crossover interference (Navarro et al. 1997; Andolfatto 
et al. 2001 ). Therefore, if two crossover events occur within 
the inversion loop created at meiosis, they are more likely 
to span the center of the inversion than segments closer to 
the inversion breakpoints (Navarro et al. 1997). Chromo- 
some arm XR differs by a single inversion between D. pseu- 
doobscura and D. persimilis, and the breakpoints of this 
inversion have been mapped (Noor et al. 2007; Bhutkar 
et al. 2008). The following markers were used to assay in- 
terspecies double-crossover rate inside the inversion — DPS 
X063 (XR_group6: 12,588,339, center of inversion), cen- 
tromeric breakpoint (XR_group6: 6,21 9,093), and telomer- 
ic breakpoint (XR_group8: 7,199,440). Primer sequences 
are included in supplementary file, Supplementary Material 
online. All individuals were genotyped at the center of the 
inversion and at least one of the breakpoint markers. Dou- 
ble crossovers were defined as individuals with an allele in 
the center of the inversion that was different from the 
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alleles at the two breakpoints, indicating exchange within 
the inverted segment. 

Calculation of Relative Node Depth between 
Non-SR D. persimilis, SR D. persimilis, and 
D. pseudoobscura 

A total of four XR-chromosome sequences were used for 
this study — one D. miranda, one D. pseudoobscura, and 
two D. persimilis. Drosophila miranda (Mt. St. Helena 22) 
and SR D. persimilis (Mt. St. Helena 39) sequences were 
generated as lllumina 76 bp paired-end reads (NCBI 
SRP006823). DNA for the SR sequence was isolated from 
a single male fly from the F 2 generation of a female 
captured in 1997 from Mount St. Helena, California, 
and lab strain male (X chromosome was wild caught), 
and amplified for genome sequencing using Qiagen's 
Whole Genome Amplification service (REPLI-g Service, 
Single Tube (1 00 |ig). Cat. no. 805923 Hilden, Germany). 
Drosophila miranda genome data was isolated from 
1 5 flies from an inbred stock. SR D. persimilis reads were 
aligned to the D. pseudoobscura reference genome v2.9 
using bwa-0.5.5 (file: dpse-all-chromosome-r2.9, align- 
ment settings -e 4 -t 7) (Li and Durbin 2009). Consensus 
assemblies for D. persimilis were generated using the 
alignments output by bwa and Samtools 0.1.6 pileup 
option (pileup settings -s -c -T 0.9 -N 1) (Li et al. 
2009). Drosophila miranda sequence reads were aligned 
using the same pileup process as SR D. persimilis and fil- 
tered through custom Python scripts that exclude indels 
with <70% of reads supporting them, 5 bp flanking pu- 
tative indels, bases with less than 3 read coverage, and 
bases with a quality score (Phred scale) less than 30. An- 
notated alignments between the released genomes of D. 
persimilis (MSH3, Clark et al. 2007) and D. pseudoobs- 
cura (MV2-25, Richards et al. 2005) were produced pre- 
viously from Kulathinal et al. (2009). The D. persimilis, SR 
D. persimilis, and D. miranda sequences were consoli- 
dated into one alignment file using D. pseudoobscura 
as a reference (available from the Dryad repository at 
doi:10.5061/dryad.7q0nq). Five megabases of both the 
centromeric and telomeric ends of this chromosome 
were excluded due to low diversity in these regions 
(Andolfatto and Wall 2003). This input file was then di- 
vided into 200-Kb windows and processed with a series 
of custom Perl scripts to both excise introns and correct 
out of frame coding sequences caused by indels in align- 
ment file. Due to genetic similarities between each of the 
sequences and D. pseudoobscura, gene annotations for 
each were assumed to be the same as the annotated 
D. pseudoobscura sequence. 

Custom Perl scripts (available from the Dryad repository 
at doi:10.5061/dryad.7q0nq) were written to calculate av- 
erage pairwise sequence difference (n), defined as the 



number of mismatches divided by the total number of bases, 
along third codon positions of nondegenerate codons (C4) 
between 200-kb windows of each of the annotated D. 
pseudoobscura and D. persimilis sequences. C4 sites were 
used because they have been shown to have less evolution- 
ary constraints than other coding regions (de Proce et al. 
2009). Relative node depth (RND; Feder et al. 2005), a mea- 
sure of species divergence that uses an outgroup to account 
for differences in mutation rate, was then calculated for 
these windows by dividing the pairwise differences of each 
sequence pair by the average of the pairwise sequence dif- 
ferences between each of the sequences and D. miranda 
(Machado et al. 2007). 

High-Throughput Genotyping to Analyze Recom- 
bination Rate Changes throughout the Genome 
Outside Inversions 

A subset of 480 individuals from the cross described in 
"Single-generation estimate of inversion crossover rate 
in hybrids" was genotyped to assay recombination rate 
changes relative to the within-species rate outside of the 
inversion due to the interchromosomal effect and at inver- 
sion boundaries (see supplementary table 1, Supplemen- 
tary Material online). This cross design did not account 
for direction of the cross because a previous study in this 
system showed no difference in ICE between Ft hybrids 
backcrossed to either D. persimilis or D. pseudoobscura 
(Ortiz-Barrientos et al. 2006). 

DNA was isolated from the 480 backcross offspring 
individually at the Genomic Sciences Lab at North Carolina 
State University for subsequent genotyping with 96 single 
nucleotide polymorphisms (SNPs) (see methods and scripts 
described in Stevison and Noor (2010)). A subset of 
42 markers corresponds to positions on the second chromo- 
some of the within-species recombination map generated 
previously for D. persimilis (Stevison and Noor 201 0). These 
markers were designed to assess the pervasiveness of the in- 
terchromosomal effect on this chromosome due to inversion 
heterozygosity in the Ft females of this cross. This chromo- 
some was targeted due to the fine-scale nature of the exist- 
ing within-species map (markers apx. 150-200 kb apart) 
relative to other chromosomes. 

Another subset of 48 markers (see supplementary table 1 , 
Supplementary Material online) was designed to further 
fine-map recombination reduction at inversion boundaries 
as in Kulathinal et al. (2009). These markers correspond 
to 2-3 Mb outside the breakpoints of each of the three ma- 
jor chromosomal inversions differentiating D. pseudoobs- 
cura and D. persimilis (six boundary regions total). Of 
these 48 markers designed for inversion boundaries, there 
were six on each XL boundary, seven on each XR boundary, 
and six on each Chromosome 2 (C2) boundary. Additionally, 
there were two markers within the XR inversion (plus two 
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additional markers designed as genotyping replicates at the 
center of the inversion, see below) and four within both the 
XL and C2 inversions. 

Finally, six markers were designed to duplicate the gen- 
otyping markers used in the inversion recombination survey 
of XR, so as to include this subset of individuals in the data 
for the original cross (see Single-Generation Estimate of 
Inversion Crossover Rate in Hybrids). 

SNP markers were screened in all offspring using the 
lllumina BeadXpress platform (Fan et al. 2003) (lllumina, 
Inc. San Diego, CA) at the Genomic Analysis Facility within 
the Duke University Center for Human Genome Variation. 
The output consisted of raw genotypes at all markers for 
all individuals. Eleven total markers were not useful for anal- 
ysis for various reasons, such as polymorphism within strains 
or monomorphism between strains. Furthermore, eight 
individuals did not amplify at any of the markers, and two 
individuals were removed from two different subsets due 
to greater than 10% missing data in that subset of markers. 
One individual was removed priorto assembly of the C2 map, 
and another individual was removed from the XL boundary 
map prior to assembly. 

The raw data was processed to assess crossovers at each 
interval via scripts from Stevison and Noor (2010). For C2, 
a Kosambi cM/Mb value was calculated and compared with 
the published rate in D. persimilis (Stevison and Noor 2010) 
to calculate a fold change difference between the hybrid 
map and the pure species map (fig. 2). Previous studies 
comparing recombination rates between D. pseudoobscura 
and D. persimilis have noted very tight correspondence 
(Ortiz-Barrientos et al. 2006; Stevison and Noor 2010); 
therefore, we limited our pure species comparison with 
the D. persimilis recombination rates. To examine how 
variation in interspecies divergence may be predicted by 
changes in recombination due to inversions, we calculated 
for each recombination interval interspecies sequence 
divergence. For consistency with our previous measures of 
divergence (see "Calculation of RND between non-SRD. per- 
similis, SRD. persimilis, and D. pseudoobscura"), we defined 
divergence as the number of mismatches between D. persi- 
milis and D. pseudoobscura divided by the total number of 
bases, along third codon positionsof nondegeneratecodons 
(C4) based on the annotated D. pseudoobscura and D. 
persimilis sequences. 



Results 

Single-Generation Estimate of Inversion Crossover 
Rate in Hybrids 

Of 9,739 individuals screened, one sample was confirmed 
as a double recombinant across the XR inversion, with the 
D. pseudoobscura allele at the center of the inversion and 
the D. persimilis allele just inside both breakpoints. Confir- 



mation consisted of first repeating the initial genotyping at 
all three markers spanning the XR inversion. Next, a gene 
conversion event was ruled out by confirming the D. pseu- 
doobscura allele at a second marker (DPSX051, chromo- 
some scaffold XR_group6: position 12953124) near the 
center of the inversion, —360 kb from the original geno- 
typing marker. This exchange is much larger than gene 
conversion tract lengths in Drosophila (Hill i ker et al. 
1994; Schaeffer and Anderson 2005). Finally, the sample 
was confirmed via Sanger sequencing at markers nearby the 
three genotyping markers, again showing a mismatch be- 
tween the allele at the center of the inversion relative to 
the inversion breakpoints. Hence, we estimate a rate of 
0.01% or ~10~ 4 double-crossover events across this 12.5 
Mb inversion in a single generation, though this figure 
may underestimate the total double-crossover rate because 
only a single position within the inversion was surveyed. 

Comparisons between Intraspecific and Interspecific 
Inversion Difference of XR Chromosomes 

We calculated Feder et al.'s (2005) RND for whole XR- 
chromosome sequences between D. persimilis (per) and 
both SR D. persimilis (SRper, comparison I: Intraspecies), 
and D. pseudoobscura (ps, comparison II: Interspecies) 
(fig. ~\A). Greater RND values indicate more differentiation 
has occurred between the focal species after correcting 
for divergence to the outgroup (D. miranda), whereas 
smaller RND values indicate less divergence between 
the focal species. Each 200-kb window was placed into 
one of four groups according to whether it was more 
(A) or less (B) than 2.5 Mb outside the XR inversion 
and more (D) or less (C) than 2.5 Mb inside the inversion 
(fig. 16). Overall, calculations of mean RND values for 
each region of comparison I resulted in the following rank 
order: A < B < D < C (fig. 1C). Hence, the regions just 
inside the inversion breakpoints had significantly higher 
relative divergence than all other regions. In contrast, 
comparison II resulted in the same rank order but with 
different regions being significantly different from one 
another: A < B < D < C (fig. 1C). Hence, the regions 
far outside the inversion breakpoints had significantly 
lower relative divergence than all other regions. As pre- 
dicted, the RND between D. persimilis — SR and D. pseu- 
doobscura (no inversion difference) was not significantly 
different for any comparison. 

Analysis of Recombination Rate Reduction in 
Single-Generation Hybrids at Inversion Boundaries 

Genotypes of markers at inversion boundaries and the 
markers within each inversion were compared for each in- 
dividual to determine how far outside the inversion break- 
points strong suppression of recombination extends (see 
detailed results in supplementary table 2, Supplementary 
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0.674 


C: 0.203 
D: 0.192 
59 
0.111 



Fig. 1. — Intraspecific versus interspecific differentiation along XR. (A) Arrows represent the orientation along XR of each of the three genomes 
compared. (6) Ideogram of the XR-chromosome arm (modified from Schaeffer et al. [2008] fig. 8D and reproduced with permission from Genetics 
Society of America) with the inverted segment highlighted by the green circle, and the segments compared labeled above the chromosome. (Q 
Summary of all pairwise comparisons of mean RND within each chromosome region (A, B, C, and D) of intraspecific (I) and interspecific (II) RND values. 
Each cell gives the mean RND for each region, total number of data points, and P value using a two-sample Wilcoxon test. 



Material online). At the C2 centromeric boundary, no re- 
combinants were recovered among 471 progeny surveyed 
as far as 2.4 Mb and as far as 2.56 Mb on the telomeric 
boundary. At the XR centromeric boundary, no recombi- 
nants were recovered as far as 2.44 Mb, and on the telomer- 
ic boundary, no recombinants were recovered as far as 2.79 
Mb. At the XL boundaries, there was an absence of 



observed recombination as far as 3 Mb on the centromeric 
side of the inversion, and no recombinants were recovered 
on the telomeric side of the inversion as far as 2.73 Mb. 

The range between minimum recombination suppression 
and the maximums presented above is summarized in table 1 
and compared with previous results from Kulathinal et al. 
(2009). The results from the current study used more 
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Table 1 

Summary of Recombination Suppression at Inversion Boundaries as Compared with Kulathinal et al. (2009) 





Kulathinal et al. (2009) (Mb) 


Current Study 


Narrowest Range (Mb) 


Centromeric side XL 


0-2.84 


2.75-3 Mb 


2.75-2.84 


Telomeric side XL 


0.4-2.8 


2.73 Mb+ 


2.73-2.8 


Centromeric side XR 


0-3.35 


2.18-2.44 Mb 


2.18-2.44 


Telomeric side XR 


2.1-2.8 


2.79 Mb+ 


2.79-2.8 


Centromeric side C2 


0-4.55 


2.32-2.56 Mb 


2.32-2.56 


Telomeric side C2 


0-2.77 


2.29-2.4 Mb 


2.29-2.4 



Note. — The refined range combines the results of the two studies. See supplementary table S2, Supplementary Materia! online for detailed results. 



markers to fine map both the lower and upper limit of this 
range. The third column represents the narrowest range as 
summarized from both studies. 

Analysis of ICE Using High-Throughput Genotyping 

The recombination rate of hybrids along C2 ranged from 
0 (within inverted regions) to 1 8.92 cM/Mb with an average 
interval size of 500 kb. The comparable recombination rate in 
the same intervals for pure species ranged from 0.1 5 to 20.54 
cM/Mb. Although the range of recombination rate values 
was similar, the distribution of these events was significantly 
different in the hybrids in 31 of 47 intervals. The ratio of cM/ 
Mb in hybrids to cM/Mb in pure species (D. persimilis) yielded 
a calculated fold change, which ranged from 0 to 9.14. The 
fold change is highest (>2-fold) in the first 4 Mb and the last 
5.5 Mb of the chromosome. Figure 2 shows a plot comparing 
recombination rates in D. persimilis and hybrids as well as 
a plot of the log (base 1 0)-normalized fold-change of recom- 
bination rate in hybrids relative to pure species. For this study, 
differences in recombination rates of hybrids along chromo- 
some 2 (bearing an inversion) can be attributed to inversion 
heterozygosity of C2 (intrachromosomal effect) and inversion 
heterozygosity of XL and XR (interchromosomal effect). 

Because previous studies on ICE have also observed 
large effects in distal and proximal regions (which nor- 
mally have low recombination in Drosophila), we tested 
for an association between recombination rate with 
D. persimilis and the fold change observed in hybrids. 
We found that regions with low recombination rate in 
D. persimilis were among those with the highest log fold 
change in hybrids (P < 0.0001 ; r 2 = 0.373; N = 36), sug- 
gesting these regions are most susceptible to changes in 
recombination. We excluded regions with no observed re- 
combination in hybrids (inside the inversion and within 3 
Mb of the inversion breakpoints). 

Next, we tested for an association of nucleotide diver- 
gence between species and hybrid recombination rate. 
We calculated divergence between D. pseudoobscura 
and D. persimilis at 4-fold degenerate third position codons 
(C4). Further, we calculated a corresponding diversity mea- 
sure along the same regions between two sequences of 
D. pseudoobscura at C4 sites, to correct for variation in lev- 
els of diversity contributing to variation in divergence. In 



a multiple regression analysis with log fold change as 
the response variable and the diversity and divergence cal- 
culations as predictor variables, we found log fold change 
explains a significant proportion of the variation in diver- 
gence between species (P = 0.0088; N = 31) but does 
not account for the variation in diversity within species 
(P = 0.0831), suggesting that the correlation observed 
was not driven by lower segregating ancestral diversity 
in chromosomal ends (Noor and Bennett 2009). In the 
above analysis, regions with no recombination in hybrids 
were excluded, along with one interval with no C4 bases 
for the diversity measure, and the last 4 intervals of the 
chromosome which are misassembled in the D. pseu- 
doobscura genome making it difficult to obtain reliable di- 
vergence estimates. When we further reduced the analysis 
to exclude chromosome ends, the results remained the 
same with a significant result for divergence (P = 
0.0049; N = 1 8), but not for diversity (P = 0.0569), despite 
the reduced variation in the log fold change variation be- 
cause ICE affects the ends most strongly. 

Discussion 

Between-Species Exchange across XR Inversion 
Detected in One Generation; However, Prolonged 
Evidence of Exchange Not Observed 

We found that the observed rate of exchange via double 
crossovers along the largest inversion (XR inversion, 12.5 
Mb) differentiating D. pseudoobscura and D. persimilis 
was similar to published rates of exchange across inversions 
segregating within species (Levine 1956; Ishii and Charles- 
worth 1977) but much lower than recombination incollinear 
regions greater than —2.5 Mb outside of the inversion or 
within the same region in homokaryotypes (Kulathinal 
et al. 2008; Stevison and Noor 2010). This observed rate 
of double-crossover in inversion heterokaryotypes ap- 
proaches the expected genomic rate of gene conversion 
(~1 0~ 5 -1 0~ 6 ), estimated using the rosy locus in Drosophila 
(Smith et al. 1 970; Chovnick 1 973). Our estimate of 1 0~ 4 is 
also much higher than estimates previously used in models 
examining the role of inversions in persistence of species with 
gene flow (Navarro etal. 1997; Kirkpatrickand Barton 2006; 
Feder and Nosil 2009). 
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Fig. 2. — Inter-chromosomal effect. (A) A comparative plot of recombination rate along the second chromosome between the published 
D. persimilis map and the map of recombination rate in between species hybrids. (6) A plot of the fold change difference in recombination rates. 



Higher Rates of Long-term Exchange within XR 
Inversion in Intraspecies versus Interspecies 
Genomic Comparisons 

One way to test if genetic background plays a role in differ- 
entiation in heterokaryotypes is to compare differentiation 
along segments of an inversion segregating within species 
with differentiation along the same segments of the same 
inversion between species. We did this by comparing the 
XR-chromosome arm of D. persimilis with D. pseudoobscura 
(heterokaryotype between species) and D. persimilis with D. 
persimilis — SR (heterokaryotype within species) and examin- 
ing how the pattern of differentiation varied between these 
comparisons. As expected, divergence relative to an out- 
group (RND) is higher in all regions in the interspecies com- 
parison between D. persimilis and D. pseudoobscura versus 
the intraspecies comparison between D. persimilis and D. 
persimilis — SR. 

The expected effect of gene flow through double crossing- 
over in inversion heterozygotes should reduce divergence at 
the midsegment of the inverted region relative to regions near 
the breakpoints (regions C vs. D in fig. 1), where double cross- 
ing-over is significantly reduced (Navarro et al. 1997). This 



prediction is consistent with our observed differences in 
RND between regions just inside the breakpoint and near 
the midsegment of the inversion, showing that differentiation 
among regions within the inverted region (C and D vs. A and B) 
persists both within- and between-species. The observation 
that the RND at region C is significantly higher than at region 
D for the intraspecies comparison (P = 0.0001 ) but not for the 
interspecies comparison (P = 0.111), suggests that higher 
rates of exchange via double crossing-over have occurred 
within species than between species. These differences are 
likely due to the higher frequency of inversion heterozygotes 
forming within species, the higher rates of exchange in the 
more closely related homologous chromosomes, and the pu- 
tative lack of selective consequences of exchange products 
within species. These results are also consistent with a recent 
population genetics analysis, which failed to find evidence for 
historical gene flow between D. persimilis and D. pseudoobs- 
cura within the same inverted segment (Stevison 201 1). 

The mean difference between regions in collinear parts 
of the chromosome and regions within 2.5 Mb outside the 
inversion (A vs. B) is significantly higher in the interspecies 
comparison (P = 0.025). However, in the intraspecies 



Genome Biol. Evol. 3:830-841. doi:10.1093/gbe/evr081 Advance Access publication August 9, 2011 



837 



Stevison et al. 



GBE 



comparison, these regions are not significantly different (P 
= 0.913), suggesting that the impact of recombination 
suppression just outside the inversion boundaries may also 
be dependent on the frequency of inversion heterozy- 
gotes. Furthermore, in comparison I, only group C showed 
a mean RND significantly different from the collinear region 
A(P< 0.0001). In comparison II, however, all other groups 
showed significantly different mean RND values from 
group A (A vs. B, P = 0.025; A vs. C, P = 0.0008; and 
A vs. D, P = 0.01 51). These combined observations suggest 
that the overall strength of heightened divergence for be- 
tween species inversions is dependent on the frequency of 
inversion heterozygotes. Finally, both the pattern of in- 
creased exchange at the midsegment relative to the break- 
point and significantly increased divergence at the 
breakpoints relative to collinear regions are consistent with 
previous studies modeling gene flow across inversions (i.e., 
Navarro et al. 1997) and experimental results (Schaeffer 
and Anderson 2005; Machado et al. 2007). 

Recombination Suppression Extends 2.5-3 Mb 
beyond Between-Species Inversion Boundaries 

One of the more puzzling phenomena associated with inver- 
sion heterozygosity is the extension of recombination sup- 
pression beyond inversions, outside the breakpoints. Here, 
we refined the known ranges of recombination suppression 
at inversion boundaries for the three major inversions differ- 
entiating D. pseudoobscura and D. persimilis (table 1). The 
results in table 1 clearly show that the level of recombination 
suppression due to inversion heterozygosity cumulatively ex- 
tends 2.5-3 Mb beyond the breakpoints of each inversion, 
adding 5-6 Mb to the total expected size of the inversion it- 
self. Because recombination is not suppressed in these bound- 
ary regions in the absence of an inversion (Stevison and Noor 
2010), it is not immediately obvious how recombination sup- 
pression would extend this far outside of the inversion loop. 
Previous studies have accredited this result to the difficulty of 
the synaptonemal complex from forming at inversion bound- 
aries (Roberts 1 976). Another possibility is that because inver- 
sion boundaries tend to accumulate in unstable/repetitive 
regions of chromosomes (Andolfatto et al. 1999; Caceres 
et al. 1999; Ranz et al. 2007), inversion boundaries could 
serve to recruit recombination suppression over the length 
of the inversion, triggered perhaps by heterozygosity imme- 
diately outside the inversion boundary. If inversion boundaries 
were indeed the molecular trigger for reduced recombination 
inside inverted regions, it would follow that these regions are 
also susceptible to recombination suppression. 

Interchromosomal Effect Is Highest in Regions of 
Low Recombination and Corresponds to Variation 
in Between-Species Divergence 

Our study investigates both intra- and interchromosomal 
effects (ICE) using —50 markers along chromosome 2. 



Because we analyze hybrids between D. pseudoobscura 
and D. persimilis, the ICE we observe is due to heterozygos- 
ity at three inversions — two on the X and one on the second 
(focal chromosome). Previous studies have observed ICE 
yielding differences in recombination rates as high as 
250% higher than standard map distance (Schultz and 
Redfield 1951), whereas we observed greater than 800% 
higher recombination rates, with an average of 224%. 

The higher proportional increase in recombination rate in 
this study relative to earlier studies could result from 1) het- 
erozygosity for more inversions (other studies observe ICE 
with only 1-2 heterozygous inversions), 2) the size of the in- 
tervals at which we analyzed recombination differences in hy- 
brids, and/or 3) higher differentiation in collinear regions due 
to our use of inversions differentiating species rather than 
populations of the same species. Previous studies averaged 
ICE over larger intervals and therefore may have masked 
the effects of smaller regions with very strong ICE. In our 
study, the average percent change in the first and last 5 
Mb is —250% higher recombination in inversion heterozy- 
gotes (more closely matching differences within species for 
single inversions); therefore, the size of intervals assayed in 
the current study is most likely responsible for the higher ob- 
served change in recombination rates. 

Several studies on inversions segregating within Dro- 
sophila melanogaster have shown that ICE has a stronger 
impact on proximal (centromeric) and distal (telomeric) 
portions of chromosomes (Schultz and Redfield 1951; 
Lucchesi and Suzuki 1968; Portin and Rantanen 1990; 
Krimbas and Powell 1992) and that these regions tend 
to have more restricted recombination within homokar- 
yotypes. Similar to previously published results, we found 
that ICE was strongest in regions distal and proximal on 
chromosome 2. Because these regions are known to have 
lower recombination rates overall, we showed that, irre- 
spective of chromosome position, regions of low recom- 
bination were most susceptible to ICE, supporting that 
these regions may be less resilient to disruptions in recom- 
bination rate. The observation of heightened effect on 
chromosome ends/centromeric regions has also been ob- 
served in maize and grasshoppers (Lucchesi and Suzuki 
1 968); however, restricted recombination in these regions 
is not universal. For example, although in Drosophila, 
both centromeric and telomeric regions tend to have re- 
duced recombination (Begun et al. 2007), in mammals, 
telomeric regions tend to have higher than average 
recombination (Kong et al. 2002). Therefore, future re- 
search could examine if ICE impacts telomeric regions 
in mammals. 

Finally, we examined how changes in recombination rates 
in hybrids correspond to patterns of gene flow and differ- 
entiation, as an extension of how inverted regions (which 
have low recombination in hybrids) often bear high diver- 
gence between species. To test for this, we calculated 
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divergence between D. pseudoobscura and D. persimilis and 
found a strong correlation with log fold change of recom- 
bination in hybrids relative to pure species. When we cor- 
rected for regions of low within species diversity (e.g., 
chromosome ends), we still observed a significant associa- 
tion, showing that low diversity at chromosomal ends 
(where ICE is strongest) does not influence the relationship 
between divergence and recombination rate changes in 
hybrids. Hence, it appears that the interchromosomal effect 
may actually "increase" interspecies gene flow outside of 
the inverted regions — a factor not considered previously 
in the effects of inversions on species persistence. 

Conclusions 

Our results suggest that models determining the role of in- 
versions in maintaining species should place less emphasis 
on single-generation recombination rates and more empha- 
sis on both the expected frequency of heterokaryotes and 
potential fitness consequences that prevent the success 
of progeny that carry exchange products from propagating 
in the next generation (see also Faria and Navarro 2010; 
Jackson 201 1). Although exchange products are rarely ob- 
served in inversion heterozygotes, the frequency of hetero- 
zygotes forming dictates the frequency of exchange 
products in the population. This frequency is a function 
of both the expected level of migration between species 
and any premating barriers to exchange. In D. pseudoobs- 
cura and D. persimilis, we know that the rate of formation of 
hybrids in sympatry is 10~ 4 (Powell 1983); however, the 
overall species hybridization rate would need to account 
for the proportion of D. pseudoobscura found in sympatry. 
Therefore, for between-species inversions, exchange is the 
product of two rare events (rare formation of inversion het- 
erozygotes and rare double-crossover events). Whereas, for 
within-species inversions, it is the product of one rare and 
one common event, with inversion heterozygote formation 
expected to be higher within-species. 

Furthermore, because the accumulation of genetic in- 
compatibilities occurs gradually over time, selection does 
not need to be very strong in early generations when migra- 
tion is nearly zero. However, upon secondary contact, higher 
rates of migration, and thus selection are likely, indicating 
that these parameters should be considered nonindepend- 
ent. These three factors — the frequency of heterokaryo- 
types/migrants, the fitness consequences of exchange in 
heterokaryotypes, and the rate at which genetic incompat- 
ibilities accumulate — are likely the most important factors 
contributing to the absence of long-term exchange de- 
tected between D. pseudoobscura and D. persimilis along 
the inverted region on the XR-chromosome arm. 

In the past 10 years, interest in the role of chromosomal 
inversions in speciation has been rekindled based on the in- 
herent properties of inversions to restrict recombination in 



heterokaryotypes. This reduced recombination protects ex- 
isting adaptive complexes and genetic incompatibilities and 
allows for the accumulation of additional incompatibilities 
between species. The research presented here confirms that 
some features found in inversions segregating within spe- 
cies apply to interspecies inversion differences but also iden- 
tifies several potential differences and hypothesizes their 
causes. Further research should explore some of the pat- 
terns suggested here, particularly considering 1) which fac- 
tors may maintain divergence in inverted regions between 
species despite detectable double crossover events, 2) the 
relationship between size of recombination intervals on 
the intensity of inter- and intrachromosomal effects (ICE) ob- 
served, and 3) how much ICE increases interspecies gene 
exchange outside inverted regions. 

Supplementary Material 

Supplementary file and table S1 and S2 are available at 
Genome Biology and Evolution online (http://www.gbe 
.oxfordjournals.org/). 

Acknowledgments 

We would like to thank S. McGaugh for preparation of sam- 
ples for genomic sequencing and alignment of sequences to 
reference genome. We thank S. Bennett and A. Boehling for 
technical support. We also thank C. Smukowski, S. 
McGaugh, and T. Mitchell-Olds for comments on the man- 
uscript. Funding for this work was provided by grants from 
the National Institutes of Health (GM076051, GM086445, 
and GM092501 ) and from the National Science Foundation 
(DEB0715484 to M.A.F.N.). 

Literature Cited 

Andolfatto P, Depaulis F, Navarro A. 2001 . Inversion polymorphisms and 
nucleotide variability in Drosophila. Genet Res. 77:1-8. 

Andolfatto P, Wall JD. 2003. Linkage disequilibrium patterns across 
a recombination gradient in African Drosophila melanogaster. 
Genetics 165:1289-1305. 

Andolfatto P, Wall JD, Kreitman M. 1999. Unusual haplotype structure 
at the proximal breakpoint of ln(2L)t in a natural population of 
Drosophila melanogaster. Genetics 153:1297-1311. 

Aquadro CF, Weaver AL, Schaeffer SW, Anderson WW. 1 991 . Molecular 
evolution of inversions in Drosophila pseudoobscura — the amylase 
gene region. Proc Natl Acad Sci USA. 88:305-309. 

Begun DJ, et al. 2007. Population genomics: whole-genome analysis of 
polymorphism and divergence in Drosophila simulans. PLoS Biol. 
5:2534-2559. 

Betran E, Rozas J, Navarro A, Barbadilla A. 1997. The estimation of the 
number and the length distribution of gene conversion tracts from 
population DNA sequence data. Genetics 146:89-99. 

Bhutkar A, et al. 2008. Chromosomal rearrangement inferred from 
comparisons of 12 Drosophila genomes. Genetics 179:1657-1680. 

Brown JD, O'Neill RJ. 2010. Chromosomes, conflict, and epigenetics: 
chromosomal speciation revisited. Annu Rev Genomics Hum Genet. 
11:291-316. 



Genome Biol. Evol. 3:830-841. doi:10.1093/gbe/evr081 Advance Access publication August 9, 2011 



839 



Stevison et al. 



GBE 



Caceres M, Ranz JM, Barbadilla A, Long M, Ruiz A. 1 999. Generation of 

a widespread Drosophila inversion by a transposable element. 

Science 285:415-418. 
Carneiro M, Blanco-Aguiar JA, Villafuerte R, Ferrand N, Nachman MW. 

2010. Speciation in the European rabbit (Oryctolagus cuniculus): 

islands of differentiation on the X chromosome and autosomes. 

Evolution 64:3443-3460. 
Carson HL. 1946. The selective elimination of inversion dicentric 

chromatids during meiosis in the eggs of Sciara impatiens. Genetics 

31:95-113. 

Chovnick A. 1973. Gene conversion and transfer of genetic information 
within inverted region of inversion heterozygotes. Genetics 
75:123-131. 

Cirulli ET, Kliman RM, Noor MAF. 2007. Fine-scale crossover rate 

heterogeneity in Drosophila pseudoobscura. J Mol Evol. 64:129-135. 
Clark AG, et al. 2007. Evolution of genes and genomes on the 

Drosophila phylogeny. Nature 450:203-218. 
Coop G, Przeworski M. 2007. An evolutionary view of human 

recombination. Nat Rev Genet. 8:23-34. 
Coyne JA, Meyers W, Crittenden AP, Sniegowski P. 1993. The fertility 

effects of pericentric inversions in Drosophila melanogaster. 

Genetics 134:487-496. 
de Proce SM, Halligan DL, Keightley PD, Charlesworth B. 2009. Patterns 

of DNA-sequence divergence between Drosophila miranda and 

D. pseudoobscura. J Mol Evol. 69:601-611. 
Dobzhansky T. 1936. Studies on hybrid sterility. II. Localization of 

sterility factors in Drosophila pseudoobscura hybrids. Genetics 

21:113-135. 

Dobzhansky T, Epling C. 1948. The suppression of crossing over in 

inversion heterozygotes of Drosophila pseudoobscura. Proc Natl 

Acad Sci USA. 34:137-141. 
Dobzhansky T, Sturtevant AH. 1938. Inversions in the chromosomes of 

Drosophila pseudoobscura. Genetics 23:28-64. 
Dobzhansky TG. 1937. Genetics and the origin of species. New York: 

Columbia University Press. 
Dobzhansky TG, Epling C. 1944. Contributions to the genetics, 

taxonomy, and ecology of Drosophila pseudoobscura and its 

relatives. Washington (DC): Carnegie Institute. 
Fan JB, et al. 2003. Highly parallel SNP genotyping. Cold Spring Harb 

Symp Quant Biol. 68:69-78. 
Faria R, Navarro A. 2010. Chromosomal speciation revisited: rearranging 

theory with pieces of evidence. Trends Ecol Evol. 25:660-669. 
Feder JL, Nosil P. 2009. Chromosomal inversions and species differences: 

when are genes affecting adaptive divergence and reproductive isolation 

expected to reside within inversions? Evolution 63:3061-3075. 
Feder JL, et al. 2005. Mayr, Dobzhansky, and Bush and the complexities 

of sympatric speciation in Rhagoletis. Proc Natl Acad Sci USA. 

102:6573-6580. 

Fitzpatrick CL, Stevison LS, Noor MAF. 2009. Fine-scale crossover rate 

and interference along the XR-chromosome arm of Drosophila 

pseudoobscura. Drosoph Inf Serv. 92:27-29. 
Foss E, Lande R, Stahl FW, Steinberg CM. 1993. Chiasma interference as 

a function of genetic-distance. Genetics 133:681-691. 
Gerton JL, et al. 2000. Global mapping of meiotic recombination 

hotspots and coldspots in the yeast Saccharomyces cerevisiae. Proc 

Natl Acad Sci USA. 97:1 1383-1 1390. 
Gruneberg H. 1935. A new inversion of the X-chromosome in 

Drosophila melanogaster. J Genet. 31:163-184. 
Hasson E, Eanes WF. 1996. Contrasting histories of three gene regions 

associated with ln(3L)Payne of Drosophila melanogaster. Genetics 

144:1565-1575. 



Hey J, Nielsen R. 2004. Multilocus methods for estimating population 
sizes, migration rates and divergence time, with applications to the 
divergence of Drosophila pseudoobscura and D. persimilis. Genetics 
167:747-760. 

Hilliker AJ, etal. 1994. Meiotic gene conversion tract length distribution within 
the rosy locus of Drosophila melanogaster. Genetics 1 37: 1 01 9-1 024. 

Hoffmann AA, Rieseberg LH. 2008. Revisiting the impact of inversions in 
evolution: from population genetic markers to drivers of adaptive 
shifts and speciation? Annu Rev Ecol Evol Syst. 39:21-42. 

Ishii K, Charlesworth B. 1977. Associations between allozyme loci and 
gene arrangements due to hitch-hiking effects of new inversions. 
Genet Res. 30:93-106. 

Jackson BC. 2011. Recombination-suppression: how many mechanisms 
for chromosomal speciation? Genetica 139:393-402. 

Joyce EF, McKim KS. 2010. Chromosome axis defects induce 
a checkpoint-mediated delay and interchromosomal effect on 
crossing over during Drosophila meiosis. PLoS Genet. 6:e1001059. 

Kirkpatrick M, Barton N. 2006. Chromosome inversions, local adapta- 
tion and speciation. Genetics 173:419-434. 

Kong A, et al. 2002. A high-resolution recombination map of the 
human genome. Nat Genet. 31:241-247. 

Krimbas CV, Powell JR. 1992. Drosophila inversion polymorphism. Boca 
Raton (FL): CRC Press. 

Kulathinal RJ, Bennett SM, Fitzpatrick CL, Noor MAF. 2008. Fine-scale 
mapping of recombination rate in Drosophila refines its correlation 
to diversity and divergence. Proc Natl Acad Sci USA. 
105:10051-10056. 

Kulathinal RJ, Stevison LS, Noor MAF. 2009. The genomics of 
speciation in Drosophila: diversity, divergence, and introgression 
estimated using low-coverage genome sequencing. PLoS Genet. 
5:e1000550. 

Laayouni H, Hasson E, Santos M, Fontdevila A. 2003. The evolutionary 
history of Drosophila buzzatii. XXXV. Inversion polymorphism and 
nucleotide variability in different regions of the second chromosome. 
Mol Biol Evol. 20:931-944. 

Leman SC, Chen Y, Stajich JE, Noor MA, Uyenoyama MK. 2005. 
Likelihoods from summary statistics: recent divergence between 
species. Genetics 171:1419-1436. 

Levine RP. 1956. Crossing over and inversions in coadapted systems. Am 
Nat. 90:41-45. 

Li H, Durbin R. 2009. Fast and accurate short read alignment with 
Burrows-Wheeler transform. Bioinformatics 25:1754-1760. 

Li H, et al. 2009. The sequence alignment/map format and SAMtools. 
Bioinformatics 25:2078-2079. 

Lucchesi JT, Suzuki DT. 1968. The interchromosomal control of 
recombination. Annu Rev Genet. 2:53-86. 

Machado CA, Haselkorn TS, Noor MAF. 2007. Evaluation of the 
genomic extent of effects of fixed inversion differences on 
intraspecific variation and interspecific gene flow in Drosophila 
pseudoobscura and D. persimilis. Genetics 175:1289-1306. 

Machado CA, Kliman RM, MarkertJA, Hey J. 2002. Inferring the history of 
speciation from multilocus DNA sequence data: the case of Drosophila 
pseudoobscura and close relatives. Mol Biol Evol. 19:472-488. 

Maynard Smith J, Maynard Smith S. 1954. Genetics and cytology of 
Drosophila subobscuraM\\\. Heterozygosity, viability and rate of 
development. J Genet. 52:152-164. 

McGaugh SE, Machado CA, Noor MAF. forthcoming. Genomic impacts 
of chromosomal inversions in parapatric Drosophila species. Trans R 
Soc B Biol Sci. 



840 



Genome Biol. Evol. 3:830-841 . doi:10.1093/gbe/evr081 Advance Access publication August 9, 201 1 



Effect of Inversions on Recombination 



GBE 



Modrich P, Lahue R. 1 996. Mismatch repair in replication fidelity, genetic 
recombination, and cancer biology. Annu Rev Biochem. 65: 
101-133. 

Navarro A, Betran E, Barbadilla A, Ruiz A. 1997. Recombination and 
gene flux caused by gene conversion and crossing over in inversion 
heterokaryotypes. Genetics 146:695-709. 

Nobrega C, Khadem M, Aguade M, Segarra C. 2008. Genetic exchange 
versus genetic differentiation in a medium-sized inversion of 
Drosophila: the A(2)/A(st) arrangements of Drosophila subobscura. 
Mol Biol Evol. 25:1534-1543. 

Noor MAF, Bennett SM. 2009. Islands of speciation or mirages in the 
desert? Examining the role of restricted recombination in maintain- 
ing species. Heredity 103:439-444. 

Noor MAF, Garfield DA, Schaeffer SW, Machado CA. 2007. Divergence 
between the Drosophila pseudoobscura and D. persimilis genome 
sequences in relation to chromosomal inversions. Genetics 177: 
1417-1428. 

Noor MAF, Grams KL, Bertucci LA, Reiland J. 2001. Chromosomal 

inversions and the reproductive isolation of species. Proc Natl Acad 

Sci USA. 98:12084-12088. 
Novitski E, Braver G. 1954. An analysis of crossing over within a heterozygous 

inversion in Drosophila melanogaster. Genetics 39:197-209. 
Ortiz-Barrientos D, Chang AS, Noor MAF. 2006. A recombinational portrait 

of the Drosophila pseudoobscura genome. Genet Res. 87:23-31. 
Pegueroles C, Ordonez V, Mestres F, Pascual M. 2010. Recombination 

and selection in the maintenance of the adaptive value of inversions. 

J Evol Biol. 23:2709-2717. 
Portin P, Rantanen M. 1990. Further-studies on the interchromosomal 

effect on crossing over in Drosophila melanogaster affecting the 

preconditions of exchange. Genetica 82:203-207. 
Powell JR. 1983. Interspecific cytoplasmic gene flow in the absence of 

nuclear gene flow — evidence from Drosophila. Proc Natl Acad Sci U 

S A. 80:492-495. 

Powell JR. 1992. Inversion polymorphisms in Drosophila pseudoobscura 
and Drosophila persimilis. In: Krimbas CV, Powell JR, editors. 
Drosophila inversion polymorphism. Boca Raton (FL): CRC Press, 
p. 73-126. 

Ranz JM, et al. 2007. Principles of genome evolution in the Drosophila 

melanogaster species group. PLoS Biol. 5:1366-1381. 
Richards S, et al. 2005. Comparative genome sequencing of Drosophila 

pseudoobscura: chromosomal, gene, and cis-element evolution. 

Genome Res. 15:1-18. 
Rieseberg LH. 2001. Chromosomal rearrangements and speciation. 

Trends Ecol Evol. 16:351-358. 
Robbins LG. 1974. Exchange within heterozygous inversions in 

Drosophila melanogaster. Genetics 77:105-114. 
Roberts PA. 1976. The genetics of chromosome aberration. In: 

Ashburner M, Novitski E, Wright TRF, editors. The genetics and 

biology of Drosophila. London: Academic Press, p. 67-184. 
Rockman MV, Kruglyak L. 2009. Recombinational landscape and 

population genomics of Caenorhabditis elegans. PLoS Genet. 

5:e1000419. 



Rozas J, Aguade M. 1 994. Gene conversion is involved in the transfer of 

genetic information between naturally-occurring inversions of 

Drosophila. Proc Natl Acad Sci USA. 91:11517-11521. 
Schaeffer SW, Anderson WW. 2005. Mechanisms of genetic exchange 

within the chromosomal inversions of Drosophila pseudoobscura. 

Genetics 171:1729-1739. 
Schaeffer SW, et al. 2008. Polytene chromosomal caps of 1 1 Drosophila 

species: the order of genomic scaffolds inferred from genetic and 

physical maps. Genetics 179:1601-1655. 
Schultz J, Redfield H. 1951. Interchromosomal effects on crossing over 

in Drosophila. Cold Spring Harb Symp Quant Biol. 16:175-197. 
Smith PD, Finnerty VG, Chovnick A. 1970. Gene conversion in 

Drosophila — non-reciprocal events at maroon-like cistron. Nature 

228:442-444. 

Spurway H, Philip U. 1952. Genetics and cytology of Drosophila 
subobscura. 7. Abnormal gene arrangements in element A. J Genet. 
51:198-215. 

Stephan W, Langley CH. 1998. DNA polymorphism in lycopersicon and 

crossing-over per physical length. Genetics 150:1585-1593. 
StevisonLS. 201 1 . Causes and consequences of recombination rate variation 

in Drosophila [dissertation]. [Durham (NC)]: Duke University. 
Stevison LS, Noor MAF. 2010. Genetic and evolutionary correlates of 

fine-scale recombination rate variation in Drosophila persimilis. J Mol 

Evol. 71:332-345. 
Strasburg JL, Scotti-Saintagne C, Scotti I, Lai Z, Rieseberg LH. 2009. 

Genomic patterns of adaptive divergence between chromosomally 

differentiated sunflower species. Mol Biol Evol. 26:1341-1355. 
Stump AD, et al. 2007. Genetic exchange in 2La inversion hetero- 
karyotypes of Anopheles gambiae. Insect Mol Biol. 16:703-709. 
Sturtevant AH. 1919. Inherited linkage variation in the second 

chromosome. In: Morgan TH, Bridges CB, and Sturtevant AH, 

editors. Contributions to the genetics of Drosophila melanogaster. 

III. Washington (DC): Carnegie Institute, p. 305-341. 
Sturtevant AH, Beadle GW. 1936. The relations of inversions in the X 

chromosome of Drosophila melanogaster to crossing over and 

disjunction. Genetics 21:554-604. 
Sturtevant AH, Dobzhansky T. 1936. Geographical distribution and 

cytology of "sex ratio" in Drosophila pseudoobscura and related 

species. Genetics 21:473-490. 
Weinstein A. 1918. Coincidence of crossing over in Drosophila 

melanogaster (Ampelophila). Genetics 3:135-172. 
Wesley CS, Eanes WF. 1994. Isolation and analysis of the breakpoint 

sequences of chromosome inversion in(3l)Payne in Drosophila 

melanogaster. Proc Natl Acad Sci USA. 91:3132-3136. 
White BJ, et al. 2009. The population genomics of trans-specific 

inversion polymorphisms in Anopheles gambiae. Genetics 

183:275-288. 

White MJD, Morley FHW. 1955. Effects of pericentric rearrange- 
ments on recombination in grasshopper chromosomes. Genetics 
40:604-619. 



Associate editor: Marta Wayne 



Genome Biol. Evol. 3:830-841. doi:10.1093/gbe/evr081 Advance Access publication August 9, 2011 



841 



